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(Dated: March 20, 2013) 

A spontaneous electric polarization in the many-body ground state of graphene and other sim- 
ilar materials, such as boron nitride and bilayer graphene, is predicted. This transition to the 
ferroelectric ground state is a consequence of strong Coulomb interaction between 7r-electrons in 
two-dimensional lattices. 



Graphene, the first truly two-dimensional (2D) mate- 
rial [lH3j], consists of a monolayer of carbon atoms ar- 
ranged in a hexagonal lattice, Fig. [TJ for a review see 
Each carbon atom has four electrons on its outer shell. 
Three of them (er-electrons) are strongly bound to the 
neighboring atoms and thus form the lattice. The forth 
electrons of each atom (7r-electrons) are relatively weakly 
bound to their hosts and can move through the whole lat- 
tice. Their motion can be described by the many-body 
Hamiltonian 
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where 



U(r) = J2Ua(r-R J ) 

3=1 
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is the periodic potential of the lattice, U a {r) is the atomic 
potential and Rj are the coordinates of the atoms. In 
the intrinsic material (which we will consider here) the 
number of 7r-electrons N equals the number of atoms. 

The Schrodinger problem with the Hamiltonian (JTJ 
is usually solved neglecting electron-electron interaction. 
Omitting the Coulomb term in ([T]) one reduces the many- 
body problem to the single-particle one for an electron in 
the periodic potential of the hexagonal lattice. The 
latter problem is solved within the tight-binding approx- 
imation Q which results in the single-particle electron 
spectrum shown in Fig. [2j The electron-electron inter- 
action is then taken into account only by assuming that 
all states of the lower energy band are occupied and all 
states of the upper band are empty. The many-body 
ground state wave function in this approach has thus 
the form of a single Slater determinant composed of the 
single-particle tight-binding orbitals. This ground state 
is non-degenerate. Similarly one treats the problem of bi- 
layer graphene or of the two-atomic monolayer hexagonal 
lattice (e.g. boron nitride). 

The goal of this Letter is to show that this commonly 
used approach is not completely valid. Consider two sorts 
of atoms, A and B, with the static atomic polarizabili- 



ties a a and Qg. Being isolated and placed in an exter- 
nal static electric field E , they acquire dipole moments 
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FIG. 1. The hexagonal two-dimensional lattice with two 
atoms per unit cell. The coordinates of the first (black cir- 
cles) and second triangular sublattice (open circles) are given 
by the vectors a(ni£i + 71.2^2) and a(ni£i + + respec- 
tively, where £1, £2, and £ are defined in ([3]) and a is the 
lattice constant. 



d,A = o>aE and <Lb — (XbEq, respectively. Now assume 
that these atoms are arranged in the plane 2D hexagonal 
lattice, Fig. [1] so that the atoms A occupy the sites of 
the triangular sublattice Ai = a(mj£i +^£2) an d the 
atoms B - of the sublattice Bi = Ai + b, b = a£, where 
a is the lattice constant, rrij and n$ are integer, and the 
vectors £1, and £2 and £ are defined as follows 



/ 1/2 \ /-1/2 
Ci = V3/2 , C2 = V3/2 

' (3) 

Consider the response of such a lattice to the external 
electric field E parallel to the layer. Under the action of 
this field each atom gets polarized but in the lattice one 
should take into account that the local field acting on the 
atom at the given point differs from the external one by 
the field produced by all other dipoles at this point @: 




FIG. 2. The tight-binding spectrum of electrons in a two- 
dimensional hexagonal lattice of Fig. [TJwith identical atoms 
in both sublattices. 
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Since all atoms in each sublattice are identical, i.e. 
df' B — d A ' B , the system of equations (HI)-© can be 
written in the form: 



d A d B 



(6) 



f d A d B } 

d B = a B < E + — ^S B a + -^tS bb > , (7) 
^ 2a i 2a 6 J 

where the sums Sa a = S bb correspond to the summa- 
tion over the dipoles of the same sublattice, 



Saa = E TX 
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and the sums Sab = Sba ~~ to the dipoles of the other 
sublattice, 
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The series in © and ([9]) converge, and their numerical 
calculation gives 



Saa ~ 11.034, S A b ~ 23.151. 



(10) 



If the atoms A and -B are identical (the case of 
graphene) then a a — ct B — a, and Eqs. © and ([7]) 
give the dipole moment of each atom 

aE 



dA = d B = 



1 - /3(Saa + Sab) 
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and the macroscopic polarizability of the lattice (per unit 
area) 



X2D,hex 



1 - p(S AA + Sab) ' 
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Here (3 = a/2a 3 and N s = 4/V30 2 is the surface den- 
sity of atoms (the polarizability \2D of a 2D lattice has 
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the dimensionality of length) . One sees that if the lattice 
constant a is large, a 3 ^> 0!(Saa + Sab)/2, the correction 
due to the local field effects is small. It is however posi- 
tive, i.e. the local field, really acting on each of the lattice 
atoms, is larger than the external field En and points in 
the same direction. If a decreases, the dipole moment d 
and the lattice polarizability X2D,hex grows and tends to 
infinity when a 3 -> a(S A A + S AB )/2. The diver gency of 
the polarizability X2D.hex physically means the instabil- 
ity of the lattice: even without the external field En the 
lattice gets spontaneously polarized since a small fluc- 
tuation of the dipole moment of any atom leads to the 
polarization of the whole system in the same direction. 
The 2D hexagonal lattice of atoms is thus stable with 
respect to the spontaneous electric polarization only if 
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Now we can check whether the stability condition (TIB")) 
is satisfied for graphene. The lattice constant in graphene 
is a = 2.46 A [4|. The atomic polarizability of the carbon 
atom is ac ~ 1.7 A 3 0. Substituting these data in (fi~3|) 
we get for graphene the number, 



ac (Q 



Sab) « 1-95 > I, 
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which is substantially above the stability boundary. The 
hexagonal graphene lattice is thus unstable with respect 
to the transition to a ferroelectric ground state with a 
finite in-plane electric polarization. The direction of the 
spontaneous dipole moment in the plane of the 2D lat- 
tice is arbitrary, so that the many-body ground state of 
graphene is degenerate, in contrast to the single-particle 
picture outlined above. 

The same conclusion should be drawn for boron nitride 
and probably for many other 2D crystals. For the two- 
atomic hexagonal lattice Eqs. (JB]) - © give the following 
expressions for the dipole moments 
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The two-atomic hexagonal lattice will then be stable if 
the denominator in Eqs. (| 1 5 1) — (jT^j) is positive. Using 
parameters of boron (as = 3.04 A 3 ) and nitrogen (ajv = 
1.10 A 3 ) 0, and the lattice constant of BN, a = 2.52 A, 
we get however a negative value of the denominator, 

(1 - p B S AA )(l - PnSaa) - PbPnS 2 ab « -1.78, (17) 

i.e. the BN lattice is also unstable and is in the ferro- 
electric ground state. One can show that the lattice of 
bilayer graphene is unstable, too. 




FIG. 3. The local field effects in (left) three and (right) two 
dimensional simple cubic/square lattice. 



The much stronger tendency of 2D lattices (as com- 
pared to 3D ones) to the ferroelectric transition has a 
simple explanation. The increase of the local field due to 
dipole-dipole interaction is also the case in 3D crystals 
and leads to the well known Clausius-Mossotti formula 
for the effective macroscopic 3D polarizability [f| 
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where N v is the volume density of atoms. In three di- 
mensions, however, this effect is not so strong. This is 
illustrated in Fig. [3] by the example of the simple cu- 
bic/square lattices. In the 3D case, Fig. [3] (left), each 
dipole is surrounded by four nearest neighbors in the az- 
imuthal direction and two neighbors in the vertical di- 
rection. Each azimuthal dipole creates the field —d/r 3 
opposite to the external one; the "north" and "south" 
dipoles produce the field +2e£/r 3 in the same direction 
as the external field. The resulting force from the near- 
est neighbors is thus zero. In the 2D case, Fig. [3] (right), 
there are only two azimuthal dipoles so that the result- 
ing field from the nearest neighbors equals +2 d/r 3 and 
enhances the external one. 

The stability condition (x > 0) is therefore typically 
not violated in 3D crystals. For example, for graphite 
parameters the denominator in (fTHJ) is positive (= 0.19). 
Whether any 2D (or ID) lattice is unstable with respect 
to the ferroelectric transition depends on the polarizabil- 
ity of constituent atoms a, the lattice period a, as well 
as on the type of the lattice. The hexagonal lattice is 
unstable, as follows from (TTBl and (fT0|) , if 



(19) 



> 0.0585. 

a 6 

The similar condition for a square 2D lattice, 
a 



> 0.2214, 



a-' 



(20) 



shows that it is more stable. In a linear (ID) chain of 
atoms a spontaneous electric polarization along the chain 
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arises if 

> 0.2080. (21) 

a 6 

We emphasize that in all cases the system gets polarized 
due to the displacement of the electron density with re- 
spect to the atomic lattice which is assumed to be fixed. 

The predicted ferroelectric transition is a consequence 
of electron-electron interaction, which is seen to be so 
strong in 2D lattices, that it leads to the reconstruction 
of the ground state (this effect could not be obtained 
by treating the Coulomb interaction as perturbation). 
The tight-binding theory of the single-particle spectrum 
of graphene electrons, Fig. [2j should thus be reconsid- 
ered. Notice that the spectrum of Fig. [2] was first calcu- 
lated by Wallace [j| who studied the properties of three- 
dimensional graphite, in which the spontaneous polariza- 
tion is not the case and hence the use of the tight binding 
theory is beyond question. 

The interaction induced spontaneous electric polariza- 
tion and the resulting degeneracy of the ground state of 
graphene and similar materials may help to better under- 
stand such effects as the minimal electric conductivity of 
graphene 0, HJ and the recently observed [§] spectrum 
reconstruction in bilayer graphene. A detailed study of 
all consequences from the predicted instability would be 
very interesting and will be done in future publications. 
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